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Abstract 

An overdamped system with a linear restoring force and two multiplicative colored 
noises is considered. Noise amplitudes depend on the system state x as x and \x\". 
An exactly soluble model of a system is constructed due to consideration of a specific 
relation between noises. Exact expressions for the time-dependent univariate proba- 
bility distribution function and the fractional moments are derived. Their long-time 
asymptotic behavior is investigated analytically. It is shown that anomalous diffu- 
sion and stochastic localization of particles, not subjected to a restoring force, can 
occur. 
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1 Introduction 



A time evolution of a stochastic system can be described by a Langevin equa- 
tion [1]. In such approach, the influence of a fluctuating environment is ac- 
counted for by means of an external noise term. It is necessary to deter- 
mine exact statistical characteristics of a system in terms of given statisti- 
cal characteristics of a noise. Assumption of Gaussian white (delta-correlated) 
noise considerably simplifies the problem [2]. The time evolution of the system 
is Markovian process [3], and its univariate probability distribution function 
(PDF) and transition probability density satisfy the Fokker-Planck equation 
[2,4,5]. In specific cases one can determine an exact time-dependent solution 
of such an equation, and a stationary PDF can be obtained in general. 

However, Gaussian white noise has some unphysical properties [6], therefore 
sometimes its application is not justified. Gaussian colored noise with an arbi- 
trary correlation function is more realistic model of a fluctuating environment. 
Although the time evolution of the stochastic system is non-Markovian pro- 
cess, in speciflc cases its univariate PDF satisfies the time dependent Fokker- 
Planck equation [7,8,9,12]. Exact statistical characteristics of undamped [10] 
and damped free particles [11] driven by additive Gaussian colored noise of 
overdamped particles in a quadratic potential driven by multiplicative Gaus- 
sian colored noise [12] are found. In particular, it is shown that anomalous 
diffusion [10,11,12] and stochastic locahzation [11,12] of free particles can oc- 
cur. 

In a more general fiuctuating environment of a system is modeled 

by two noise sources. Most commonly, cross-correlated Gaussian white noises 
are considered. The appropriate Fokker-Planck equation is derived and the 
two-noise Langevin equation can be reduced to a stochastically equivalent 
one-noise Langevin equation [13,14,15]. However, models which demonstrate, 
for example, nonequilibrium fiuctuation-induced transport [16,17,18], double 
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stochastic resonance [19], require to consider a colored noise source and a 
white noise source. The Ornstein-Uhlenbeck process with exponential cor- 
relation function is widely used as Gaussian colored noise. But the corre- 
sponding Fokker-Planck equation is derived using approximate methods, for 
instance, conventional small-r theory [20], the decoupling theory [21], the uni- 
fied colored noise approximation [22] , etc. Within this approximations, systems 
driven by two Ornstein-Uhlenbeck processes are studied [23,24], the noise self- 
correlations and the cross-correlation have the same correlation time. There 
are no general methods for obtaining exact statistical characteristics of sys- 
tems with two colored noises. Therefore, specific exactly soluble models are 
needed. 

In this paper, we generalize a nonlinear system [12] with colored noise to the 
case of two colored noises and study the possibihty of an exactly solvable 
model constructing. The system state parameter x{t) evolves according to the 
following Langevin equation: 

x{t) + [K + h{t)]x{t) = \x{t)rf2{t), (1) 



where x{0) — xo{> 0), k{> 0) and a{< 1) are real-valued parameters, fi{t) 
and f2{t) are colored noises with zero mean and known statistical character- 
istics. The assumption of positivity of x{0) does not restrict the generality of 
the model (1). 

We note that for < a < 1 the solution of Eq. (1) is not unique at x = 0. In 
this paper, the point a; = is considered to be a regular point and the solution 
of Eq. (1) coincides with the solution of the equation 

{x{t) + [k + Mt)]x{t)}\x{t)r = f2{t) [x{Q) =xo>0] (2) 



for all times t > 0. Although one can obtain an exact solution of Eq. (2) 
for the case of arbitrary independent noises fi{t) and /2(t), it is not possible 
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to determine exact statistical characteristics of x{t). Therefore, we consider a 
specific relation between noises fi{t) and f2{t), which allows to obtain exactly 
the time-dependent univariate PDF and the fractional moments of x{t). We 
find analytically their long-time asymptotics and show that in the special case 
K — anomalous diffusion and stochastic localization of particles can occur. 

The paper is organized as follows. In Sec. II, we solve exactly Eq. (2) with ar- 
bitrary independent noises fi{t), f2{t) and consider a specific relation between 
noises. In the same section we obtain the time-dependent univariate PDF and 
the fractional moments of x{t). Their long-time asymptotic behavior is studied 
analytically in Sec. III. Our results are summarized in Sec. IV. 



2 General analysis 



To construct an exactly soluble model of a system described by the Langevin 
equation (2), we obtain its exact solution for the case of arbitrary independent 
noises fi(t) and /2(i)- By introducing the new variable y{t) — x{t)\x{t)\~°' , the 
nonlinear stochastic differential equation (2) is transformed to the linear one 

y{t) + [a; + (1 - = (1 - a)Mt), (3) 



where u; — {1 — a)K and y{0) — xl The solution of Eq. (3) is given by (see, 
for example, Ref. [25]) 



y{t) = e 



-ujt-{l-a)Fi{t) 



4"" + {l-a) J (iT/2(r)e""+(^-")^i(") 



(4) 



where 



F,{t) = J drMr). 



(5) 
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Substituting Eq. (4) into expression x{t) — y{t)\y{t)\'^/^^ °'\ we can write the 
solution of Eq. (2) as 



where 



xl-'^ + il-a) I dTf2{r)e' 



u}T+{l-a)Fi{T) 



(7) 



To obtain the time-dependent univariate PDF of x{t), we should determine 
the univariate PDFs of Fi{t) and F{t). The univariate PDF of Fi{t) is easily 
found if is Gaussian noise with zero mean and the arbitrary correlation 
function 

{h{t)Mt'))^n{\t-t'\). (8) 



Indeed, by virtue of the linear dependence Fi{t) on fi{t), Fi{t) is Gaussian 
process with zero mean and the dispersion 

t 

al {t)^2 j dun {u){t-u). (9) 



The simplest way to obtain the univariate PDF of F{t) is to consider the 
specific relation between noises fi{t) and f2{t)- We assume that 

m = /(i)e(--^)^^W, (10) 



where f{t) is Gaussian noise with zero mean and the arbitrary correlation 
function 

U{t)f{t'))=r{\t-t'\). (11) 



5 



Then F{t) is Gaussian noise with the mean m{t) and the dispersion a'^{t) 
m{t) = 4""e-'"*, (12) 

-cot * 

a^{t) = 2(1 - a) / dur{u)smh[u}{t - u)]. (13) 

K J 


Thus, the Langevin equation (2) with arbitrary noises /i(t) and f2{t) is re- 
duced to the foUowing Langevin equation: 

{x{t) + [K + Mt)]x{t)}\x{t)r = /(t)e("-^)^^W [x{0) =Xo> 0], (14) 

where fi{t) and f{t) are statistically independent Gaussian colored noises with 
known statistical characteristics. 

The random process p{t) — F(i)|F(i)|"/(^~") is power-normal [12]. Let Pp{p, t) 
denote the univariate PDF that p{t) = p, 

'^^^''^-V2^a(t)bh^"n 2a^W T ^ ^ 

The random process l{t) = e^^*^*^ is logarithmic- normal [26]. Let Pi{l,t) denote 
the univariate PDF that l{t) — I, 



In^Z 



2al{t)\ 



(16) 



for Z > and Pi{l,t) — for I < 0. According to Eq. (6), the solution of 
Eq. (14) can be rewritten as x{t) = p{t)/l{t), where p{t) and l{t) are statis- 
tically independent power-normal and lognormal processes. It is not difficult 
to get the conclusion, that by introducing the new variable p{t) — x{t)l{t), 
Eq. (14) can be obtained from the Langevin equation with one noise 

m+f^pimpit)r-fit) [p(o)=po>o], (i?) 



where po — xq. 
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Let Px{x, t) be the univariate PDF that x{t) — x. It can be calculated as (see, 
for example, Ref. [27]) 



oo 

Px{x,t) — J duuPi{u,t)Pp{ux,t). 



(18) 



Substituting Eq. (15) and Eq. (16) into Eq. (18), we finally obtain the time- 
dependent univariate PDF of x{t) 

1 — a /■ , f In^M [ux\x\~" — m{t)]'^ ^ 



Px{x,t) = 



27iaa{t)a{t)\x\ 



du exp 



2^^(t) 



2a\t) 



,(19) 



where Uait) = {1 — a)ai{t). 

We note if /i(t) = [o"^(t) for all times t] the small neighborhood of the 
point u = 1 gives basic contribution and Eq. (19) is reduced to the PDF of 
the power-normal distribution (15). If f{t) = [cr^{t) for all times t] we 
take into account 



lim 



a{t) 



\x 



l-a 



X 



a{t) 



\x 



l-a 



exp < — 



u 



1 



u — 



m(t) 



x\x\ 



m{t) 



x\x\ 



(20) 



where 5{x) is the Dirac delta function. Eq. (19) is reduced to the PDF of the 
lognormal distribution [12] 



Px{x,t) 



27rai{t)x 



exp 



2af{t) V xo 



In \- Kt 



(21) 



for X >0 and Px{x,t) = for a; < 0. Taking into account Eqs. (20) and (21), 
we note that Eq. (19) reduces to the delta function S{x — xq) at time t — 0. 

We can now determine the time-dependent fractional moments of x{t). They 
are defined as 



oo 

m;:(t) = J dxPx{x,t)\x\''-"x\ 



(22) 
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where r is a real number, and v = or 1. We substitute Eq. (18) into Eq. (22). 
By introducing the new variable y — ux, one can write 



(23) 



We take into account that the second integral is the well-known fractional 
moments of power-normal distribution [12]. Calculating the first integral, we 
finally obtain the fractional moments of x{t) 



m'^it) = ^^a^-\t) exp \r^(rl{t)/2 - a\t)/4 



x{D.^[-a{t)] + {-irD.Mt)]}. 



(24) 



Here ^ = 1 + r / (1 ~ a) , a{t) = m{t)/a{t), and D^^{z) is the integral represen- 
tation of the Weber parabolic cylinder function [28] 



^2/4 °^ 



D_^{z) = ^ / dyy^-'e 



(e > 0), 



(25) 



where r(^) is the gamma function. If ^ < all fractional moments diverge. 
Since mQ{t) — 1, Px{x,t) is properly normalized. 

We note the asymptotic of 'm'^{t) as (7f{t) — > for all times t is obtained 
directly from Eq. (24) 

Kit) = ^e-«^W/V^-i(t){i?_,[-«W] + {-irD_,[am- (26) 



Eq. (26) is the fractional moments of the power- normal distribution [12]. In 
order to determine the fractional moments as cr^{t) — > for all times t, we use 
the asymptotic form 



lim (7^ ^(i)exp 



4a2(t) 



m(t) 



ait) 
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oo 



m 



m 




5{s±l). 



(27) 







Substituting Eq. (27) into Eq. (24), we obtain 



'm^{t) — XqBxp -r^al{t) — rut . 



(28) 



Eq. (28) is the fractional moments of the logarithmic- normal distribution [12]. 
3 Asymptotic behavior 

In the previous section, we have constructed an exactly soluble model with two 
multiplicative Gaussian colored noises. We have obtained exact expressions 
for the time-dependent univariate probability distribution function and the 
fractional moments. In this section, we determine the long-time behavior of 
particles, i — > oo. Cases k > and k — Q are considered separately. 

3.1 K > 

In this case m(oo) = and a(oo) = 0. The asymptotic form of Eq. (13) at 
i = oo is given by 



Since r{u) — > as — > oo, (7^(oo) < oo. The asymptotic form of Eq. (9) can 
be written as 



oo 




(29) 







al{t) ~ 2tRi{t) {t 



oo). 



(30) 



9 



where Ri{t) — /q ditri(ii). According to Ref. [11], if 
Ri{t) = o{l/t) {t oo), 



(31) 



then (Tj{oo) < oo, and if < i?i(oo) < oo or i?i(oo) = [-Ri(oo) is the noise 
intensity], but Eq. (31) does not hold, then cr^(oo) = oo. Thus, there are two 
quahtatively different cases of the long-time asymptotic behavior of Eq. (9), 
namely, erf (oo) < oo and crl{oo) — oo. 

To determine the asymptotic behavior of m^(i) as i ^ oo we use the formula 
^-.(0) = 2«/^-^^. (32) 



which follows from Eq. (25). Substituting Eq. (32) into Eq. (24), we obtain 

m;(oo) = £ipl2«/V«-^(oo)e^'-?(-)/2i±tl)! (^ > q). (33) 
\/2tt 2 



If Eq. (31) is fulfilled [crf{oo) < oo], all fractional moments (33) with r > a — 1 
are finite, and the stationary PDF Pst{x) = Px{x, oo) is given by 



Pstix) 



a 



27r(7a(oo)(7(oo)|a; 



oo 



du exp 



21 |2(l-a) 



U \X 



In^ju) 

'2al(oo) 2(72(oo) 



(34) 



Note that Pstix) is even, as is also reflected by m^(oo) = 0. The value of Pstix) 
at the point x = depends on a, if < a < 1 then Pst(O) = oo, if a = then 
P^^(O) = e'^?(~)/V[V27f(7(oo)], and if a < then P,t(0) = 0. 

If Eq. (31) is not fulfilled [erf (oo) = oo], then m°(oo) = oo and m^(oo) = 0. 
The stationary PDF Pgtix) — Pxix, oo) (19) can be written as 



Pstix) 



27r(7(oo)|x|°-^2ii?i(i) 



X j duexjp 



In'^iu) 



^2|^|2(l-a) 



4:{l-a)HRiit) 2(72(oo) 



it oo). (35) 
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For x 7^ 0, the first term of the exponent in Eq. (35) can be neglected compared 
to the second one, and we find 

1 



Pstix) 



L^n = hm 



i\x\^TTtRi{t) 



0. 



(36) 



For X — 0, Eq. (35) is reduced to 
PstiO) 



(l-n)exp[(l-n)^f/?,(f)] 

hm ;= — : — — ■ = oo. 

t oo 
X -» 



27r(7(oo)|x|° 



(37) 



To gain more insight into the behavior of Pgtix) in the vicinity of x 
define the probabihty 



W,{t) = J dxP^{x,t), 



0, we 



(38) 



that x{t) e (—£,£). Substituting Eq. (35) into Eq. (38), we obtain 
1 



oo 




erf 




j dz exp 

-oo 


2tA{t)_ 


V2a((X))_ 



(i^oo), (39) 



where A{t) — 2(1 — ayRi{t). To evaluate ^^(oo), we introduce the new 
variable z — sy/i, and take into account that 



lim erf 

t— >oo 



\/2(j(oo) 



h{s), 



(40) 



where h{s) is the Heaviside function. Substituting Eq. (40) into Eq. (39), we 
obtain We{oo) — 1/2. Thus, in the case (7i(oo) = oo no well-defined stationary 
PDF exists. 



3.2 K = Q 



In this subsection, we study the statistical behavior of free particles (particles 
not subjected to a restoring force [12]) under infiuence of two colored noises. 
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As mentioned in Section 1, free particles driven by additive colored noise 
or multiplicative Gaussian colored noise can display anomalous diffusion and 
stochastic localization. 



In this case (k = 0) m{t) — xl for all times t, and Eq. (13) is reduced to 

t 

a\t) = 2(1 - a f J dur{u){t - u). (41) 



Eq. (41) is qualitatively similar to Eq. (9), and its long-time asymptotic has 
the form 

a'^(t) 2(1- a)HR(t) (t^ oo), (42) 

where R{t) = j'^dur{u). If 

R{t) = o{l/t) {t oo), (43) 

then (7^(oo) < oo, and if < R{oo) < oo or i?(oo) — [R{oo) is the 
noise intensity], but Eq. (43) does not hold, then (7^(oo) = oo. So, there are 
two different cases of the long-time asymptotic behavior of Eq. (41), namely, 
cr^(oo) < oo and cr^(oo) = oo. 

It is not difficult to verify that if cr^ (oo) = oo and (T^(oo) < oo or (T^(oo) = oo 
then m°(oo) = oo and m^(oo) = 0, P^(0, oo) = oo and Px{x, oo)\^_^q = 0, 
We(oo) = 1/2. So, in this cases no well-defined stationary PDF exists. 

If (T^(oo) < oo and cr^ (oo) < oo, all fractional moments (24) with r > a — 1 
are finite 



m'^it) = ^a^~\oo) exp [rV?(oo)/2 - a\oo)/4: 



x{D_^[-a(oo)] + (-l)^D_^[a(oo)]}, (44) 
a(oo) = xJ"°/cr(oo), and the stationary PDF is given by 



12 



oo 



ln^{u) [ux\x 




} 

(45) 



Pst{x) 



1 — a 



I 



du exp 



—a 



27r(Ti(oo)(7(oo)|j; 



a 







2<7^(oo) 



2(72(oo) 



In contrast to the case k > 0, where the stationary PDF (34) is even, the 
stationary PDF (45) is not even, as is also reflected by the fractional moments, 
m^(oo) 7^ 0. 

Thus, if Eqs. (31) and (43) arc fulfiUed, then Gaussian colored noises f\{t) 
and f{t) lead to stochastic localization of free particles. This phenomenon 
was first described for free particles driven by additive [11] and multiplicative 
[12] colored noise. 

If (7^(oo) = oo and crf{oo) < oo, then m°(oo) = oo and mj(oo) = 0, 
Px{x,oo) — and W£(oo) = 0. Consequently, the stationary PDF does not 
exist, Gaussian colored noise f{t) gives rise to diffusive behavior of free par- 
ticles. This case is qualitatively similar to one considered in Ref. [12]. It 
is not difficult to verify, if < R{oo) < oo, then cr^(i) oc i as i — > oo, 
m°(i) oc and ml{t) oc t^/"^"^. The dispersion of the particle position, 

'^x(^) — (^^(^)) ~ can be rewritten as (jl{t) — m2{t) — [m\{t)Y. Hence, 

(j1{t) oc as i ^ oo, and the conditions a = 0, a < 0, and < a < 1 

correspond to normal diffusion, subdiffusion, and super diffusion, respectively. 

4 Conclusions 

We have constructed an exactly solvable nonlinear model of an overdamped 
system with two multiplicative Gaussian colored noises and studied its statis- 
tical properties. Starting from the exact solution of the Langevin equation for 
the case of arbitrary independent noises, we have considered the specific rela- 
tion between noises. In such approach, the time evolution of the system state 
has been represented as the ratio of the power-normal process to the lognor- 



13 



mal one, and the two-noise Langevin equation can be reduced to the one-noise 
Langevin equation. We have obtained exact expressions for the time-dependent 
univariate probabihty distribution function of the system state and the frac- 
tional moments. Analyzing their analytically found long-time asymptotics, we 
have shown that anomalous diffusion and stochastic localization of particles, 
not subjected to a restoring force, can occur. 
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